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I.  InttoductLoa 

A.  General: 

Energetic  materials  are  substances  containing  metastable  molecules  enable  of  undergoing 
very  rapid  highly  exothermic  reactions  [Iyer  and  Slagg,  1988].  They  are  generally  divided  into 
three  major  classes:  explosives,  propellants  and  pyrotechnics  [Meyer,  1987;  Iyer  and  Slagg, 

1988].  Our  focus  in  this  chapter  will  be  primarily  on  the  first  two  categories  (explosives  and 
propellants). 

The  study  of  energetic  systems  by  theoretical  methods  has  accelerated  dramatically  over  the 
course  of  the  last  decade  and  has  provided  considerable  insight  into  the  understanding  of  factors 
affecting  their  behavior.  In  this  chapter  we  will  discuss  first-principles  calculations  that  have  been 
carried  out  for  energetic  molecules  of  a  variety  of  structural  types,  including:  a)  nitroaromatics;  b) 
aliphatic  and  alicyclic  nitramines;  c)  highly  strained  cagelike  systems;  and  d)  high-nitrogen-content 
molecules.  Our  scope  will  include  ab  initio  molecular  orbital  methods  and  also  density  functional 
procedures.  A  specific  emphasis  will  be  the  understanding  of  factors  influencing  shock/impact 
sensitivities  and  overall  stabilities. 

Before  launching  our  discussions  of  calculations  carried  out  for  the  different  structural 
types,  it  is  useful  to  point  out  a  molecular  feature  that  is  common  to  all  of  them,  namely  a  relatively 
high  degree  of  charge  separation.  For  example,  a  component  of  many  currently-used  and 
proposed  energetic  molecules  (e.g.  1  -  6)  is  the  NO2  group  [Meyer,  1987;  Iyer  and  Slagg,  1988; 
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Bulusu,  1990].  Its  presence  introduces  significant  charge  separation,  or  local  polarity,  even  in  a 
molecule  which  has  a  zero  dipole  moment,  due  to  the  nitro  group's  strongly  electron  attracting 


2 


tendencies  [Exner,  1988].  Thus,  molecules  such  as  1  -  6  can  be  viewed  as  possessing  relatively 
high  degrees  of  local  polarity.  We  will  discuss  the  significance  of  this  for  energetic  molecules  in 
later  sections. 


B.  Structure-Sensitivity  Relationships: 

A  continuing  major  theme  of  energetic  materials  research  has  been  the  effort  to  better 
understand  the  relationship  between  the  molecular  structures  of  energetic  compounds  and  their 
sensitivities  to  specific  stimuli,  such  as  impact  and  shock.  Indeed,  this  has  been  the  driving  force 
for  many  computational  studies.  A  key  long  range  goal  is  the  design  and  manufacture  of  less 
sensitive,  less  vulnerable  munitions. 

Sensitivity  has  been  defined  as  the  ease  with  which  a  substance  subjected  to  external 
stimuli,  such  as  shock,  impact  or  heat,  can  be  made  to  undergo  detonation  [Iyer  and  Slagg,  1988]. 
An  ideal  explosive  would  be  one  which  achieves  high  performance,  but  is  insensitive  enough  to 
permit  safe  handling;  however  it  should  not  be  so  insensitive  that  detonation  becomes  difficult  or 
impractical. 

Many  kinds  of  sensitivity  have  been  identified  in  terms  of  the  nature  of  the  stimulus  causing 
detonation  [Iyer  and  Slagg,  1988];  some  of  the%  include  heat,  fiiction,  impact,  shock  and 
electrostatic  charge.  Two  of  the  most  commonly  used  measures  of  an  explosive's  tendency  to 
detonate  are  impact  and  shock  sensitivity.  Impact  tests  (of  which  there  are  several,  similar  in 
principle  but  differing  in  detail)  involve  subjecting  a  sample  to  the  intact  of  a  standard  weight 
failing  from  varying  heights  [Kamlet  and  Adolph,  1979];  the  sensitivity  is  inversely  proportional  to 
the  height  at  which  there  is  a  50%  probability  of  causing  an  explosion  (called  impact  drop  height  or 
hso).  Shock  sensitivity  is  often  measured  by  a  standard  small  scale  gap  test  [Storm,  Stine  and 
Kramer,  1990],  and  is  interpreted  as  being  directly  proportional  to  the  maximum  gap  width 
through  which  a  given  shock  wave  generated  by  another  explosive  (e.g.  RDX)  will  detonate  the 
compound  50%  of  the  time. 

Kamlet  and  Adolph's  (1979)  stoichiometric  means  of  assessing  impact  sensitivity  through 
correlations  with  "oxygen  balance"  has  been  widely  applied  in  the  area  of  energetic  materials 
research.  The  term  "oxygen  balance",  or  OBioo,  refers  to  the  amount  of  oxygen  available  for  the 
combustion  of  the  C  and  H  atoms  in  a  molecule  to  CO  and  H2O,  and  is  defined  for  C-H-N-0 
compounds  in  eq.  (1). 


OB  100  = 


100  (2no  -  nH  -  2nc  -  2ncoo) 
M 


(1) 


no,  nn  and  nc  represent  the  numbers  of  atoms  of  the  respective  elements  in  the  molecule,  ncoo  is 
the  number  of  carboxyl  groups,  and  M  is  the  molecular  weight.  Kamlet  and  Adolph  demonstrated 
that  for  families  of  high  energy  molecules  with  similiar  decomposition  mechanisms  (e.g.  trinitro- 


methyl,  N-nitro  and  gem-dinitro  compounds  and  nitroaromatics  with  and  without  a-C-H  linkages) 
there  are  £^proximately  linear  relationships  between  log  hso  and  OBioo-  In  these  correlations,  log 
h5o  increases  as  OBioo  decreases  for  any  one  family  of  molecules;  the  more  impact-sensitive 
explosives  are  those  with  the  more  positive  OBioo  values.  Stine  [1989]  has  since  proposed  a 
sensitivity  index  based  on  the  geometric  properties  of  a  regular  tetrahedron  which  has  been  shown 
to  correlate  with  OBioo  values  [Storm,  Stine  and  Kramer,  1989]. 

Although  such  stoichiometric  relationships  have  proven  to  be  quite  useful  for  a  number  of 
classes  of  compounds,  taken  separately,  they  do  not  allow  one  to  distinguish  between  structural 
isomers,  which  can  have  quite  different  sensitivities.  This  consideration  has  prompted  a  variety  of 
investigations  aimed  at  elucidating  the  relationship  of  molecular  structure  to  impact/shock 
sensitivity.  Considerable  attention  has  been  given  to  the  idea  of  a  key  "trigger"  linkage  in 
explosives,  as  was  proposed  by  Kamlet  and  Adolph  [1981].  In  light  of  numerous  studies,  both 
theoretical  and  experimental,  supporting  the  importance  of  C-NO2  and  N-NO2  bonds  in 
determining  the  sensitivities  of  many  nitroaromatics,  nitroaliphatics  and  nitramines,  it  is  the  general 
consensus  that  in  many  such  molecules,  C-NC)2  or  N-NO2  bonds  are  the  actual  "triggers". 

Table  1  lists  the  sensitivities  for  some  common  explosives.  For  the  limited  number  of 
compounds  included  in  the  table,  it  can  be  seen  that  the  nitrate  esters  1  and  2  and  the  nitramines  4 
(RDX)  and  7  (HMX)  are  more  sensitive  to  shock  and  impact  than  the  nitroaromatics  (3, 5  and  8  - 
10).  The  sensitivities  of  the  trinitroaromatics  cover  a  relatively  wide  range  and  appear  to  be  highly 
dependent  upon  the  substituents  on  the  aromatic  ring.  In  this  chapter,  factors  affecting  the  impact 
and/or  shock  sensitivities  of  explosives  will  be  examined  separately  for  different  types  of 
molecules,  when  possible,  and  later  from  a  broader  perspective,  encompassing  all  of  them. 

NO2  O2N 


7  8  9  10 

n.  Methodologies 
A.  General: 

The  calculations  to  be  discussed  in  this  chapter  involve  primarily  ab  initio  molecular  orbital 
and  density  functional  procedures.  Theoretical  background  and  computational  details  pertaining  to 
these  two  first-principles  approaches  are  reviewed  elsewhere  [ab  initio:  Hehre,  Radom,  Schleyer 
and  Pople,  1986;  Hinchliffe,  1988;  density  functional:  Parr  and  Yang,  1989;  Kryachko  and 
Ludena,  1990;  Labanowski  and  Andzelm,  1991].  The  emphasis  in  this  chapter  is  upon  the  results 
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of  these  computations.  There  have  also  been  semi-empirical  studies  of  energetic  molecules  (e.g. 

Alster,  Slagg,  Dewar,  Ritchie  and  Wells,  1981;  Turner  and  Davis,  1984),  and  some  of  these  will 
be  mentioned. 

B.  Structures: 

A  complete  theoretical  determination  of  the  equilibrium  structure  of  a  molecule  involves 
minimization  of  its  energy  with  respect  to  each  independent  geometrical  parameter  [Hehre,  Radom, 
Schleyer  and  Pople,  1986],  resulting  in  an  "optimized"  structure.  The  resulting  geometry  is 
necessarily  somewhat  dependent  upon  the  methodology  and  basis  set  used  to  perform  the 
optimization. 

Table  2  compares  some  calculated  ab  initio  Hartree-Fock  (HF),  post-Hartree-Fock  and 
density  functional  structural  data  with  experimental  values  (when  available)  for  several  small 
molecules  which  can  be  viewed  as  building  blocks  of  larger  explosives.  Several  significant  trends 
emerge  which  are  important  in  selecting  computational  procedures  and  basis  sets  for  geometry 
optimizations  of  molecules  containing  C-NO2  and  N-NC)2  linkages.  (The  term  "basis  set"  refers  to 
the  functions  used  to  build  up  the  molecular  orbitals.  The  ones  that  will  be  mentioned  are 
designated  as  STO-3G,  3-21G  and  6-31G*;  these  increase  in  size  and  flexibility  in  the  order 
given.)  Thus,  although  HF/STO-3G  structures  are  often  quite  satisfactory.  Table  2  shows  that 
N-0  bond  lengths  obtained  at  this  computational  level  are  too  large  by  about  0.05  A.  The 
HF/3>21G  N-0  distances  are  much  better,  being  overestimated  by  ^proximately  0.02  A,  while  the 
HF/6-31G*  tend  to  be  underestimated  by  about  tiw  same  amount.  The  inclusion  of  electron 
correlation  via  post-Hartree-Fock  or  density  functional  procedures  is  seen  to  yield  N-0  bond 
lengths  that  are  generally  comparable  to  the  HF/3-21G. 

Thus,  Table  2  shows  that  satisfactory  nitro  group  geometries  can  be  computed  both  with 
and  without  the  inclusion  of  electron  correlation.  In  the  firameworir  of  ab  initio  HF  theory,  it  is 
advisable  to  carry  out  geometry  optimizations  with  the  more  flexible  basis  sets  (e.g.  3-2 IG  or 
6-3 IG*).  The  effect  of  computational  method  and  level  upon  the  structural  features  of  high- 
nitrogen-content  molecules  will  be  discussed  in  detail  in  a  later  section. 

C.  Energies: 

Calculated  energies  can  be  used  to  determine  relative  stabilities  of  structural  isomers, 
reaction  energetics  and  isodesmic  reaction  energies.  Isodesmic  reactions  provide  insight  into 
anomalous  energy  effects  in  molecules  [Hehre,  Ditchfield,  Radom  and  Pople,  1970;  Hehre, 

Radom,  Schleyer  and  Pople,  1986].  They  are  hypothetical  chemical  processes  in  which  the 
number  of  bonds  of  each  formal  type  remains  the  same  on  each  side  of  the  equation,  but  their 
mutual  relationships  are  changed.  For  example,  an  isodesmic  reaction  involving  cyclopropane 
(C3H6)  is:  3C2H<5  — >  C3H6  +  3CH4.  The  AE  values  for  reactions  such  as  this  are  interpreted  as 


representing  any  deviations  from  bond  energy  additivity  and  therefore  indicative  of  any  anomalous 
energetic  effects  associated  with  the  molecule  being  investigated,  e.g.  strain,  resonance 
stabilization,  steric  effects,  etc.  (AE  <  0  indicates  the  presence  of  stabilizing  factors,  while  AE  >  0 
implies  destabilization.) 


D.  Reaction  Energetics: 

The  elucidation  of  the  decomposition  mechanisms  of  energetic  molecules  is  of  considerable 
importance,  both  for  understanding  how  the  stored  energy  is  released  and  also  for  seeking  to 
improve  the  stabilities  of  such  systems.  Theoretical  investigations  of  reaction  energetics 
complement  experimental  studies  in  this  area.  Several  of  computational  studies  will  be  discussed, 
in  which  the  interest  was  primarily  in  obtaining  reaction  energies  and/or  activation  barriers,  and 
which  were  carried  out  with  either  ab  initio  or  density  functional  procedures. 


E.  Molecular  Properties: 

1.  The  Electrostatic  Potential,  V(r): 

The  electrostatic  potential  is  a  real  physical  property  that  reflects  the  net  electrical  effect  of 
the  nuclei  and  electrons  of  a  molecule  at  any  point  in  space  r;  it  is  defmed  rigorously  by  eq.  (2). 


A 


J  Ir’-ii 


(2) 


Za  is  the  charge  on  nucleus  A,  located  at  Ra,  and  p(r)  is  the  electronic  density  function,  which  we 
determine  computationally.  The  electrostatic  potential  is  well-established  as  an  effective  guide  to 
chemical  reactivity  [Scrocco  and  Tomasi,  1978;  Politzer  and  Daiker,  1981;  Politzer  and  Truhlar, 
1981;  Politzer  and  Murray,  1991].  In  those  regions  where  V(r)  is  negative,  the  effects  of  the 
electrons  dominate  and  these  are  accordingly  attractive  sites  for  the  initial  attack  of  electrophiles. 
Predicting  sites  reactive  toward  nucleophiles  is  less  straightforward  due  to  the  effects  of  the  nuclei, 
but  we  have  shown  that  when  V(r)  is  plotted  on  the  molecular  surface,  the  regions  of  most  positive 
potential  do  allow  predictions  for  nucleophilic  processes  [Sjoberg  and  Politzer,  1990;  Politzer  and 
Murray,  1991;  Murray  and  Politzer,  1991]. 


2.  Atomic  Charges  and  Vrnid' 

We  have  recently  modified  and  extended  an  earlier  rel^onship  between  impact  sensitivity 
of  nitroaromatics  and  Vmid.  the  electrostatic  potential  at  the  midpoint  of  the  C-NO2  bonds  [Owens, 
Jayasuriya,  Abrahmsen  and  Politzer,  1985;  Murray,  Lane,  Politzer  and  Bolduc,  1990] ;  the  latter  is 
approximated  by 
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V  H  -  ^ 

'  0.5  R  ^  0.5  R 


(3) 


Qc  and  Qn  are  the  atomic  charges  on  the  carbon  and  nitrogen,  obtained  from  the  Mulliken 
population  analysis  of  the  SCF-MO  wave  function  [Mulliken,  1955];  R  is  the  C-N  bond  length  in 
bohrs.  The  relationship  between  nitroaromatic  in^ru^t/shock  sensitivities  and  VnUd  will  be 
discussed  in  a  later  section. 


3.  The  Average  Local  Ionization  Energy,  l(r): 

We  have  also  introduced  another  calculated  property  that  is  an  effective  tool  for  analyzing 
certain  aspects  of  chemical  reactivity.  This  is  the  average  local  ionization  energy,  i(r),  defmed  by 
eq.  (4)  [Sjoberg,  Murray,  Brinck  and  Politzer,  1990]: 


i(r)  =  X 

\ 


PjCr)  Ujl 
P(r) 


(4) 


Pi(r)  is  the  electronic  density  of  the  i*  molecular  orbital  at  the  point  r,  e;  is  the  orbital  energy,  and 
p(r)  is  the  total  electronic  density.  i(r)  can  be  interpreted  as  the  average  amount  of  energy  required 
to  remove  an  electron  at  any  point  r  in  the  space  of  the  molecule.  We  have  shown  that  i(r)  plotted 
on  molecular  surfaces  is  a  very  good  indicator  of  sites  for  electrophilic  attack;  diese  are  where  i(r) 
has  its  lowest  values,  indicating  the  most  easily  removed  and  hence  most  reactive  electrons 
[Sjoberg,  Murray,  Brinck  and  Politzer,  1990;  Brinck,  Murray  and  Politzer,  1991].  For  example, 
for  a  series  of  monosubstituted  benzenes,  the  surface  i(r)  minima  (Is, min)  provide  a  quantitative 
measure  of  the  activating/deactivating  and  the  directing  tendencies  of  the  various  substituents,  and 
also  correlate  well  with  the  Hammett  constants  and/or  Cm  [Sjoberg,  Murray,  Brinck  and 
Politzer,  1990].  Surface  !(r)  can  also  be  used  successfully  to  predict  the  strengths  of  acids 
[Brinck,  Murray,  Politzer  and  Carter,  1991;  Brinck,  Murray  and  Politzer,  1991].  For  24  acids  of 
various  types  (carboxylic,  inorganic,  etc.),  there  is  an  excellent  linear  relationship  between  the 
experimentally-determined  pKa  of  each  acid  and  the  Is, min  of  its  conjugate  base.  It  is  particularly 
noteworthy  that  these  relationships  are  between  a  measured  solution  property  and  a  gaseous  phase 
calculated  one. 


m.  Nitroaromatics 

Due  to  the  importance  of  nitroaromatics  as  a  class  of  explosives,  there  has  been 
considerable  interest  in  their  properties,  mechanisms  of  decomposition,  and  factors  affecting  their 
sensitivities.  In  this  section  we  will  begin  with  a  discussion  of  structural  and  electronic  properties 
of  nitroaromatics,  followed  by  a  review  of  calculations  relating  to  the  decomposition  of  some 
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members  of  this  class.  Finally  we  will  survey  relationships  found  between  calculated  properties 
and  the  impact/shock  sensitivities  of  tlMse  compounds. 

Certainly  a  relevant  feature  of  aromatics  in  general  and  one  to  which  we  will  refer  again  is 
the  influence  that  any  particular  substituent  or  group  of  substituents  can  have  upon  the  aromatic 
ring.  This  is  presumably  due  to  the  presence  of  the  polarizable  n  electrons  and  the  possibility  of 
conjugation,  and  is  reflected  in  the  fact  that  both  resonance  and  inductive  effects  play  prominent 
roles  in  determining  the  distributions  of  charge  in  these  systems.  In  this  context,  it  is  noteworthy 
that  the  role  of  resonance  is  determined  by  difference  [Exner,  1988],  often  using  a  general 
relationship  of  the  type  Cp  -  oi  =  <Tr,  where  the  Hammett  constant  Op  is  defined  in  terms  of  para- 
substituted  aromatics  (e.g.  benzoic  acids),  for  which  resonance  interactions  by  resonance  donors 
can  be  quite  significant,  and  Oi  is  determined  independently  for  series  such  as  substituted  acetic 
acids  for  which  inductive  effects  predominate. 

Nitrobenzene  (11)  can  be  described  in  terms  of  resonance  structures  such  as  11a  -  lie.  A 
number  of  theoretical  studies  have  analyzed  the  inductive  and  resonance  (i.e.  11c  -  lie) 


11  11a  11b  11c  lid  lie 


contributions  to  the  strongly  electron-withdrawing  tendency  of  the  NC)2  [Lipkowitz,  1982;  Hiberty 
and  Ohanessian,  1984;  Politzer,  Lane,  Jayasuriya  and  Domelsmith,  1987;  Ritchie,  1988].  The 
general  consensus  is  that  there  is  only  a  small  degree  of  conjugation  between  the  nitro  group  and  the 
aromatic  ring  in  nitrobenzene.  This  is  consistent  with  the  well-established  inductive  and  resonance 
substituent  constants  CTi  and  Or  of  the  nitro  group  (0.67  and  0.10,  respectively  [Exner,  1988]),  and 
with  the  relatively  low  rotational  barrier  around  the  C-NC)2  bond  in  nitrobenzene,  around  3  kcal/ 
mole  [H0g,  Nygaard  and  Sorensen,  1971;  Correl,  Larsen  and  Pedersen,  1980].  Even  in  the  case  of 
the  observed  "through-resonance"  interactions  in  some  para-substituted  nitrobenzenes  [Hammett, 
1970],  such  as  para-nitroaniline,  12,  there  is  experimental  [Lipkowitz,  1982]  and  theoretical 
[Hiberty  and  Ohanessian,  1984]  evidence  that  12c  is  the  dominant  resonance  contributor. 
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The  strongly  electron-withdrawing  nitro  group  significantly  deactivates  the  ring  an  aromatic 
ring  toward  electrophihc  aromatic  substitution  [Allinger,  Cava,  DeJongh,  Johnson.  Lebel  and 
Stevens,  1971].  This  is  reflected  in  the  electronic  density  distributions  of  nitroaromatics  and  also 
in  related  properties  such  as  the  electrostatic  potential  V(r)  and  the  average  local  ionization  energy 
!(r).  For  example,  the  regions  of  negative  electrostatic  potential  that  are  characteristically  found 
above  and  below  the  planes  of  the  benzene  ring  (Figure  1),  and  are  attributed  to  the  7C  electrons,  are 
completely  eliminated  upon  the  substitution  of  a  single  nitro  group  (Figure  2)  [Politzer,  Abrahmsen 
and  Sjoberg,  1984]. 

The  deactivating  mera-directing  influence  of  NO2  upon  electrophilic  aromatic  substitution 
can  be  predicted  from  the  calculated  average  local  ionization  energy  i(r)  of  nitrobenzene  plotted  on 
the  molecular  surface  [Sjoberg,  Murray,  Brinck  and  Politzer,  1990;  Murray,  Brinck  and  Politzer, 
1992].  The  locations  of  the  Is, min  above  the  meta  positions  indicate  that  these  are  the  positions 
where  the  most  reactive  electrons  are  found,  and  thus  are  the  most  favorable  sites  for  electrophilic 
attack.  The  magnitudes  of  these  Is, min  are  greater  than  those  of  benzene,  consistent  with  the 
deactivating  influence  of  the  NO2.  As  mentioned  earlier,  we  have  shown  that  the  Is  jnin  above  the 
para  and/or  meta  positions  of  a  series  of  substituted  benzenes  correlate  well  with  the  respective 
Hammett  constants  Op  and/or  Cm-  This  permits  us  to  predict  Hammett  constants  for  substituents 
for  which  these  are  not  known,  e.g.  the  NF2  and  N(N02)2  groups  [Sjoberg,  Murray,  Brinck  and 
Politzer,  1990;  Politzer,  Murray,  Seminario  and  Miller,  1992]. 

The  commonly  used  nitroaromatic  explosives  are  benzene  derivatives  which  contain  three  or 
more  nitro  groups  (e.g  see  5  - 10  in  Table  1).  The  sensitivities  of  these  compounds  cover  a  wide 
range  [Kamlet  and  Adolph,  1979;  Storm,  Stine  and  Kramer,  1990],  and  have  been  shown  to  be 
highly  dependent  upon  the  natures  of  the  additional  substituents  upon  the  aromatic  ring  [Owens, 
1984;  Owens,  Jayasuriya,  Abrahmsen  and  Politzer,  1985;  Murray,  Lane,  Politzer  and  Bolduc, 
1990].  Except  for  the  sensitizing  effect  of  the  OH  group  [Kamlet  and  Adolph,  1979],  which  we 
will  discuss  subsequently,  electron-donating  substituents  (e.g.  NH2,  OCH3  and  CH3)  are  generally 
desensitizing;  this  is  exemplified  by  the  aminotrinitrobenzenes  (e.g.  see  9  and  10  in  Table  1).  In 
contrast,  electron-withdrawing  groups  (including  another  NO2)  increase  sensitivity.  Owens  [1984], 


for  exanq>le,  has  demonstrated  for  a  group  of  eleven  trinitroaromatics  that  impact  sensitivity 
correlates  inversely  with  the  total  donating  ability  of  the  non-nitro  substituents  on  the  ring. 

One  interpretation  of  these  observations  is  that  resonance  donors  strengthen  to  some  degree 
the  C-NO2  bonds,  as  is  shown  below  for  2,4,6-trinitroaniline  (9): 
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As  we  have  pointed  out  earlier,  however,  the  degree  of  resonance  electron-withdrawal  by  the  nitro 
group  into  the  C-NO2  bond  region  is  believed  to  be  small,  even  with  strong  electron  donors 
present,  and  the  resonance  structures  9b  -  9d  may  play  only  a  minor  role. 

Consistent  with  an  emphasis  upon  C-NO2  bond  strength  are  the  relationships  found 
between  impact  sensitivity  and  Vmid,  the  electrostatic  potential  at  the  midpoint  of  the  C-NC)2  bond 
as  approximated  from  atomic  charges  by  eq.  (2)  [Owens,  Jayasuriya,  Abrahmsen  and  Politzer, 
1985;  Murray,  Lane,  Politzer  and  Bolduc,  1990].  In  general,  impact  sensitivity  increases  as  the 
largest  Vmid  of  the  tiinitroaromatic  molecule,  Vnud,max.  increases.  In  view  of  the  relationship 
between  eq.  (2)  and  a  bond  energy  expression  developed  by  Fliszar  (1980),  it  can  be  inferred  that 
ti»  instability  of  an  individual  C-NO2  linkage  may  be  a  key  factor  in  initiating  decomposition 
induced  by  impact.  The  presence  of  electron-donating  groups  such  as  NH2,  C)CH3  and  CH3  tends 
to  diminish  Vnud,max.  suggesting  that  these  may  reduce  the  instabilities  of  the  C-NO2  linkages 
relative  to  those  of  tiinitrobenzene. 

The  OH  group  is  an  exception  to  this  pattern.  Even  though  it  is  electron-donating, 
hydroxynitroaromatics  such  as  13  - 15  have  anomalously  high  sensitivities.  Politzer,  Seminario 


and  Bolduc  [1989]  have  investigated  the  possibility  that  this  may  be  due  to  nitronic  acid  formation. 
Using  a  model  system,  o-nitrophenol  (16),  the  proton  transfer  to  form  the  nitronic  acid  17,  shown 
in  eq.  (5),  was  predicted  (at  the  HF/3-21G  level)  to  occur  to  a  small  but  significant  extent.  It  was 
suggested  that  the  formation  of  unstable  nitronic  acids,  accompanied  by  the  disruption  of 


OH  OH  OH 


aromaticity,  may  contribute  to  the  observed  sensitizing  effects  of  hydroxyl  substituents  in 
nitroaromatics  [Politzer,  Seminario  and  Bolduc,  1989]. 

In  a  later  computational  study  [Murray,  Lane,  Politzer,  Bolduc  and  McKenney,  1990], 
hydrogen  transfer  reactions  for  two  widely-studied  nitroaromatic  explosives,  2,4,6-trinitro-toluene 
(8,  TNT)  and  l,3,5-triamino-2,4,6-trinitrobenzene  (18,  TATB),  were  investigated,  again  using 

CH3 

NO2 

18  19  20 


model  systems:  o-nitrotoluene  (19)  and  o-nitroaniline  (20).  The  objective  was  to  gain  insight  into 
the  very  different  thermal,  shock  and  impact  sensitivities  of  8  and  18  [Kamlet  and  Adolph,  1982]. 
TNT  is  thermally  labile,  self-igniting  at  temperatures  between  200  and  210°  C,  and  has 
intermediate  impact  and  shock  sensitivities.  TATB,  on  the  other  hand,  is  stable  to  temperatures 
well  above  300  °  C  and  is  so  insensitive  to  impact  and  shock  that  it  is  impractical  as  an  explosive. 
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Our  results  for  o-nitrotoluene  (19)  indicate  that  intramolecular  hydrogen  transfer  to  form  a 
bicyclic  tautomer  21  is  highly  favored  thermodynamically  at  the  HF/3-21G  level,  as  is  the  further 
loss  of  water  to  form  anthranil  22  [eq.  (6)].  An  interesting  feature  of  these  results  is  that  for  the 
model  system  19,  we  do  not  find  a  nitronic  acid  tautomer  analogous  to  17,  but  rather  the  thermo¬ 
dynamically  favored  bicyclic  systems  21  and  22.  In  this  context,  it  is  noteworthy  that  an  earlier 
semi-empirical  investigation  of  the  thermal  decomposition  of  TNT  using  1-nitropropene  (23)  as  a 
model  [Turner  and  Davis,  1984],  and  also  at  the  cJj  initio  HF/STO-3G  level  using  TNT  itself  [Cox 
and  Hillier,  1988],  found  the  formation  of  the  nitronic  acids  24  and  25  to  be  thermodynamically 
highly  unfavorable  [eqs.  (7)  and  (8)]. 


NO2  NO2 


8  25 

Evidence  of  the  formation  of  an  intramolecular,  five-membered  ring  analogous  to  22  has 
been  reported  in  low-temperature  studies  of  TNT  by  Rogers  [1967]  and  by  Dacon  et  al  [1970]. 
We  have  suggested  from  our  results  for  the  model  system  19  that  energy  produced  by  forming 
products  such  as  21  and  22  could  help  to  stimulate  further  reactions  [Murray,  Lane,  Politzer, 
Bolduc  and  McKenney,  1990].  This  idea  is  consistent  with  experimental  results  indicating  that  a 
product  of  TNT  resulting  from  a  methyl  C-H  bond  rupture  accumulates  to  a  threshold  level  and 
then  catalytically  initiates  its  continuing  exothermic  thermochemical  decomposition  [Shackelford, 
Beckmaim  and  Wilkes,  1977]. 
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For  o-nitroaniline  (20),  tautomerization  to  a  nitronic  acid  was  also  found  to  be 
thermodynamically  unfavorable  at  the  HF/3-21G  level,  while  conversion  to  benzofurazan  26  with 
the  elimination  of  water  [eq.  (9)]  was  calculated  to  be  thermochemically  nearly  neutral.  Five- 
membered  ring  formation  giving  fiirazans  (analogous  to  26)  has  been  reported  in  thermal 
decomposition  and  underwater  shock  and  impact  studies  of  TATB  (18)  [Sharma,  Hoffsomer, 
Glover,  Coffey,  Santiago,  Stolovy  and  Yasuda,  1983;  Sharma,  Forbes,  Coffey  and  Liddiard, 
1987] .  It  has  been  predicted  that  the  steps  producing  furazan  derivatives  of  TATB  [e.g.  27  in  eq. 
(10)]  are  thermodynamically  nearly  neutral,  and  this  has  been  suggested  as  helping  to  account  for 
its  very  low  sensitivity  [Storm,  Stine,  and  Bulusu,  1990].  Our  HF/3-21G  results  for  the  model 
system  20  are  consistent  with  this  view. 
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As  we  have  mentioned  earlier,  the  presence  of  the  nitro  group  in  an  organic  molecule 
introduces  some  degree  of  charge  separation,  or  local  polarity,  due  to  the  strongly  electron- 
attracting  power  of  the  NC)2;  this  occurs  even  if  the  symmetry  of  the  molecule  is  such  that  the 
overall  dipole  moment  is  zero.  The  electrostatic  potentials  shown  in  Figs.  1  and  2  for  benzene  and 
nitrobenzene  provide  some  indications  of  the  charge  separation  within  these  two  molecules. 

Indeed  V(r),  defined  by  eq.  (2),  is  a  physical  manifestation  of  the  distribution  of  electronic  and 
nuclear  charge  within  a  system. 

To  provide  a  quantitative  measure  of  the  local  polarity  that  may  be  present  even  in  a  zero- 
dipole-moment  molecule,  we  have  recently  introduced  a  polarity  index  IT,  defined  in  terms  of  the 
electrostatic  potential  V(r)  by  eq.  (11)  [Brinck,  Murray  and  Politzer,  1992]: 


S 


V(r)-V, 


dS 


(11) 


The  integral  represents  the  total  amount  by  which  the  electrostatic  potential  on  the  molecular  surface 
deviates  from  its  average  value,  Vj;  this  is  divided  by  the  surface  area  A  in  order  to  permit 
comparisons  between  molecules  of  different  sizes. 

A  computationally  practical  form  of  eq.  (1 1)  is  given  by  eq.  (12), 


(12) 


where  Vi(r)  is  the  potential  at  the  i*  point  on  the  surface  and  Vs 


1  “ 

-  y  Vi(r).  As  n  increases, 
"  i=l 


eq.  (12)  approaches  eq.  (11).  We  have  proposed  eq.  (12)  as  the  working  formula  for  11,  which 
can  be  viewed  as  the  average  deviation  of  the  electrostatic  potential  on  the  molecular  surface 
[Brinck,  Murray  and  Politzer,  1992].  We  have  computed  11  for  molecules  of  a  variety  of  types, 
and  demonstrated  its  relationship  to  an  empirical  polarity-polarizability  parameter  and  also  to  the 
dielectric  constant,  an  experimentally-determined  bulk  property. 

Table  3  lists  our  calculated  11  values  for  a  variety  of  molecules  (including  nitrobenzene,  p- 
dinitrobenzene  and  a  number  of  other  monosubstituted  aromatics),  computed  at  the  HF/STO-5G 
level  on  molecular  surfaces  defined  in  accordance  with  Bader  et  al  [1987]  by  the  0.(X)1  electrons/ 
bohr^  contour  of  the  electronic  density  [Brinck,  Murray  and  Politzer,  1992].  The  general  lack  of 
correspondence  between  11  and  the  dipole  moment  (i  can  be  seen  by  inspection  of  Table  3,  as  can 
also  the  general  relationship  between  11  and  the  dielectric  constant  E.  For  the  nitro  derivatives,  a 
general  trend  is  that  11  increases  as  the  carbon  to  NO2  ratio  decreases;  the  11  values  for  C6H5NC)2i 
/7-C6H5(N02)2  attd  CH3NO2  are  12.3, 16.5  and  19.9  kcal/mole,  respectively.  For  the 
monosubstituted  benzenes,  the  11  values  range  from  4.6  to  12.3  kcal/mole,  with  the  largest  value 
being  that  of  nitrobenzene. 

These  observations  and  the  known  general  trends  in  the  impact/shock  sensitivities  of 
nitroaromatics  led  us  to  investigate  the  possibility  that  there  might  be  a  relationship  between  11  and 
impact  sensitivity  [Murray,  Lane,  Brinck  and  Politzer,  1992].  In  Table  4  are  listed  experimentally- 
determined  impact  sensitivities  and  our  earlier  reported  Vnjicynax  values  for  a  group  of  nitro¬ 
aromatics,  along  with  their  11  values  and  surface  potential  maxima.  For  these  molecules,  excluding 
the  monohydroxy  derivative  (picric  acid,  13),  there  is  a  linear  correlation  coefficient  of  0.92 
between  impact  drop  height  and  Vnud,inax-  the  same  group,  the  impact  sensitivity  and  n  show 
a  correlation  coefficient  of  0.98  [Murray,  Lane,  Brinck  and  Politzer,  1992]. 

An  interesting  feature  of  the  surface  potentials  of  the  nitroaromatics  listed  in  Table  4  is  that 
the  most  positive  regions  of  electrostatic  potential  on  the  molecular  surface  are  generally  located 
above  the  C-NCh  bonds  and  the  center  of  the  aromatic  ring.  (This  can  be  seen  in  Fig.  3  for  1,3,5- 
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trinitrobenzene.)  These  positive  regions  are  indicative  of  sites  favorable  for  the  initial  approach  of 
nucleophiles,  and  are  consistent  with  the  observed  reactions  of  nitroaromatics  with  hydroxide  and 
alkoxide  ions  to  form  Meisenheimer  complexes  [Gibson  and  Crampton,  1979;  Crampton  and 
Wilson,  1980;  Bacaloglu,  Bunton  and  Cerichelli,  1987],  as  well  as  their  reactions  with  hydrogen 
halides  [Nielsen,  Chafln  and  Christian,  1984].  The  monohydroxy  derivative,  13,  stands  out  in 
that  the  positive  surface  potential  associated  with  its  hydroxyl  hydrogen,  about  58  kcal/mole,  is 
fully  19  kcal/mole  greater  than  that  above  its  ring  center,  and  is  in  fact  much  larger  than  the 
potential  maxima  normally  found  near  alcohol  and  phenol  hydrogens,  which  range  from  25  to  49 
kcal/mole  [Murray  and  Politzer,  1992].  The  highly  positive  nature  of  this  hydrogen  is  consistent 
with  our  earlier  study  suggesting  that  nitronic  acid  formation  through  intramolecular  proton  transfer 
may  contribute  to  the  observed  sensitizing  effect  of  OH  groups  in  nitroaromatics  [Politzer, 
Seminario  and  Bolduc,  1989]. 

rv.  Nitramines 

Energetic  secondary  nitramines  are  an  important  class  of  explosives,  as  is  seen  from  the 
fact  that  two  of  the  most  effective  (yet  highly  sensitive)  military  explosives  are  nitramines:  1,3,5- 
trinitro-l,3,5-triazacyclohexane  (4,  RDX)  and  l,3,5,7-tetranitro-l,3,5,7-tetraazacyclooctane  (7, 
HMX).  In  this  section  we  wiU  begin  with  a  discussion  of  structural  and  electronic  properties  that 
are  unique  to  nitramines,  followed  by  an  overview  of  calculations  of  decomposition  reactions  of 
some  specific  nitramines  and  perspectives  on  the  sensitivities  of  members  of  this  class  of 
explosives. 

Organic  amines  and  aza  nitrogens  are  characterized  by  the  presence  of  a  "lone  pair"  of 
electrons  which  give  rise  to  strongly  negative  regions  in  the  molecular  electrostatic  potentials; 
however  these  are  significantly  weakened  by  the  substitution  of  electron-withdrawing  nitro  groups 
on  the  amine  nitrogens  [Politzer  and  Murray,  1991].  As  an  extreme  example,  the  negative  potential 
in  the  lone  pair  region  of  dimethylamine  (28)  reaches  a  value  of  -100  kcal/mole  (at  the  HF/STO-5G 
computational  level)  [Murray,  Redfem,  Seminario  and  Politzer,  1990],  whereas  in  dimethylnitra- 
mine  (29),  there  is  qq  negative  potential  associated  with  the  amine  nitrogen  [Murray,  Redfem, 
Lane,  Politzer  and  Wilier,  1990]. 
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An  interesting  feature  of  aza  nitrogens  is  an  anomalous  stabilization  that  is  associated  with 
their  presence,  and  which  we  attribute  to  a-conjugation  of  their  lone  pairs  [Murray,  Redfem, 
Seminario  and  Politzer,  1990].  We  have  identified  this  effect  by  means  of  the  well-established 


isodesmic  reaction  procedure  [Hehre,  Ditchfield,  Radom  and  Pople,  1970;  Hehre,  Radom, 
Schleyer  and  Pople,  1986]. 

Secondary  nitramines  contain  the  subunit  R1R2N-NO2,  the  simplest  example  being  di- 
methylnitraminc,  29;  some  others  are  shown  below  (30  -  36).  In  the  course  of  several  comput¬ 
ational  studies  of  nitramines  [Politzer,  Sukumar,  Jayasuriya  and  Ranganathan,  1988;  Murray, 
Sukumar,  Ranganathan  and  Politzer,  1990;  Murray,  Redfem,  Lane,  Politzer  and  Wilier,  1990; 
Habibollahzadeh,  Murray,  Redfem  and  Politzer.  1991],  consistent  trends  have  emerged.  For 
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example,  a  decrease  in  the  degree  of  pyramidal  character  of  the  amine  nitrogen  parallels  a  decrease 
in  N-NO2  bond  length.  This  observation  is  consistent  with  greater  delocalization  of  the  amine  lone 
pair  into  the  bond  with  increasing  planarity,  as  shown  in  eq.  (13).  For  instance,  the  calculated 
HF/3-21G  N-NO2  bond  lengths  in  1,3-dinitroimidazolidine  (33)  are  0.04  A  shorter  than  those  in 
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the  less  planar  1,4-dinitropiperazine  (34)  [Murray,  Redfem,  Lane,  Politzer  and  Wilier,  1990].  In 
accordance  with  these  observations,  we  have  shown  that  the  N-NO2  rotational  barriers  in  nitra¬ 
mines  are  generally  greater  than  C-NC)2  rotational  barriers  [Habibollahzadeh,  Murray,  Grice  and 
Politzer,  1992],  and  that  the  former  are  largely  dependent  upon  the  degree  of  pyramidal  character 
of  the  amine  nitrogen  in  the  ground  state  molecule  [Habibollahzadeh,  Murray,  Redfem  and 
Politzer,  1991]. 

Due  at  least  in  part  to  their  importance  as  energetic  materials,  there  is  a  great  deal  of  interest 
in  the  stabilities  and  the  thermal  decomposition  mechanisms  of  nitramines.  The  overaU  consensus 


is  that  rupture  of  the  N-NO2  bond  is  a  key  step  in  the  process  of  decomposition,  altlK)ugh  there 
may  be  important  competing  pathways  which  may  even  predominate  in  some  instances  [Brill  and 
Oyumi,  1986;  Zhao,  Hintsa  and  Lee.  1987;  Stewart  and  Jefhies,  Zellweger,  McMillan  and 
Golden,  1989;  Saxon  and  Yoshimine,  1989;  Grodzicki,  Seminario  and  Politzer,  1990a;  Grodzicki, 
Seminario  and  Politzer,  199(%);  Habibollahzadeh,  Grodzicki,  Seminario  and  Politzer,  1991]. 

In  several  recent  computational  studies,  we  have  focused  upon  some  of  the  possible  initial 
steps  in  the  unimolecular  thermal  decompositions  of  1.3-diazacyclobutane  (37),  1,3-dinitro- 
diazacyclobutane  (38),  1,3,5-triazacyclohexane  (39)  and  its  nitramine  derivative  RDX  (4) 
[Grodzicki,  Seminario  and  Politzer,  1990;  Seminario,  Grodzicki  and  Politzer,  1991;  Grodzicki, 
Seminario  and  Politzer,  1991;  Habibollahzadeh,  Grodzicki,  Seminario  and  Politzer,  1991].  In  the 
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first  of  these,  a  vibrational  analysis  carried  out  at  the  ab  initio  HF/3-21G  level  pointed  to  two 
modes  that  clearly  support  bond  rupture  and  molecular  dissociation  [Grodzicki,  Seminario  and 
Politzer,  1990].  One  is  a  ring  mode,  while  the  other  involves  a  stretching  of  the  N-NO2  bonds;  the 
corresponding  bond-breaking  processes  are  shown  in  eqs.  (14)  and  (15): 
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An  important  difference  between  the  reactions  depicted  in  eqs.  (14)  and  (15)  is  that  the  former  has 
an  activation  barrier  while  the  latter,  N-NO2  bond  scission,  does  not.  The  energetics  analysis  in 
this  study  led  us  to  the  conclusion  that  ring  fragmentation  is  roughly  competitive  with  N-NO2  bond 
breaking  as  an  initial  step  in  the  thermal  decomposition  of  38. 

We  have  also  investigated  the  nitro-nitrite  rearrange^  nent  of  38,  given  in  eq.  (16),  and  have 
found  it  to  be  less  likely  than  the  processes  shown  in  eqs.  (14)  and  (15)  [Grodzicki,  Seminario  and 
Politzer,  1990].  Rearrangements  m,aiogous  to  eq.  (16)  are  believed  to  play  important  roles  in  the 
thermal  decoiiq)ositions  of  dimethylnitramine,  29  [Stewart,  Jewries,  Zellweger,  McMillen  and 


Golden,  1989]  and  nitramide  (H2NNO2)  [Saxon  and  Yoshimine,  1989;  Seminario  and  Politzer, 
1992]. 
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The  pathway  depicted  in  eq.  (14)  is  analogous  to  the  concerted  symmetric  ring  opening  that 
has  been  proposed  as  the  "dominant  primary  charmel"  in  the  thermal  decoti^sihon  of  RDX,  4 
[Zhao,  Hintsa  and  Lee,  1987].  We  have  recently  investigated,  using  ab  initio  and  local  density 
functional  (LDF)  methods,  the  activation  barriers  to  the  concerted  gas-phase  ring  fragmentations  of 
l,3,S-triazacyclohexane  (39)  and  its  trinitro  derivative  4,  shown  in  eqs.  (17)  and  (18) 
[Habibollahzadeh,  Grodzicki,  Seminario  and  Politzer,  1991].  We  found  that  the  calculated  barriers 
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decrease  as  the  initio  computational  level  increases,  with  the  LDF  values  being  somewhat 
smaller  in  magnitude  than  the  best  Hartree-Fock  or  post-HF.  This  conclusion  was  also  reached  in 
a  coiiq)arison  of  ab  initio  and  LDF  activation  energies  for  the  process  shown  in  eq.  (14) 
[Grodzicki,  Seminario  and  Politzer,  1991],  as  well  as  for  some  small  molecule  isomerizations 
[Seminario,  Grodzicki  and  Politzer,  1990].  Our  predicted  activation  energies  for  eqs.  (17)  and 
(18)  are  in  the  range  72  -  75  kcal/mole.  These  results  are  consistent  with  the  experimental  findings 
by  23iao  et  al  [1987]  in  molecular  beam  work  that  the  activation  energy  for  the  reaction  given  in  eq. 
(18)  is  not  more  than  80  kcal/mole. 


Understanding  factors  affecting  the  sensitivities  of  nitramines  has  been  a  key  emphasis  in 
energetic  materials  research.  As  we  have  mentioned,  considerable  attention  has  focused  upon 
cleavage  of  the  N-NO2  bond  as  a  key  step  in  the  process  of  decomposition,  initiated  by  heat,  shock 
or  impact.  We  have  recently  demonstrated  [Politzer,  Murray,  Lane,  Sjoberg  and  Adolph,  1991] 
th^  a  good  correlation  can  be  obtained  between  the  shock  sensitivities  of  the  nitramines  4, 7, 33  - 
35, 40  and  41  and  molecular  properties,  by  (a)  explicitly  taking  each  N-NO2  bond  into  account, 
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and  (b)  recognizing  that  its  effectiveness  in  triggering  decomposition  is  "diluted"  as  the  molecule 
becomes  larger.  A  straightforward  formulation  of  these  ideas  is. 

Sensitivity  .  ^  =  X  J  R,  =  (19) 

i=l  K.  i&I 

in  which  n  is  the  number  of  N-NO2  linkages  in  the  molecule,  having  N-N  distances  Ri,  and  M  is 
the  molecular  weight  M  is  taken  as  a  measure  of  molecular  size,  and  the  reciprocal  of  R,  as  our 
measure  of  bond  strength.  For  the  molecules  listed  above,  a  linear  relationship  exists  between 
(n/M)  Rave  and  shock  sensitivity,  with  a  correlation  coefficient  of  0.94.  (For  a  group  of  five 
nitroaliphatics,  a  similar  relationship  was  found,  where  n  is  the  number  of  C-NO2  linkages  and  Ri 
is  the  C-N  bond  length.) 

It  should  be  mentioned  that  for  this  group  of  nitramines,  shock  sensitivity  also  correlates 
well  with  n/M  alone;  indeed  the  latter  has  been  investigated  as  an  indicator  of  impact  sensitivity 
[Adolph,  Holden  and  Cichra,  1980].  However  the  extensive  experimental  evidence  pointing  to 
N-NO2  scission  as  a  key  step  in  nitramine  decomposition  argues  that  the  strengths  of  these  bonds 
should  be  explicitly  reflected  in  any  measure  of  sensitivity. 

Another,  yet  consistent  formulation  of  these  concepts  is  our  measure  of  local  polarity,  11 
[Brinck,  Murray  and  Politzer,  1992] .  Analogously  to  the  C-NO2  systems  in  Tables  3  and  4,  IT 
increases  as  the  ratio  of  N-NC)2  groups  to  carbon  atoms  in  the  firamework  increases,  and  so  does 
sensitivity  to  shock.  Our  calculated  IT  values  for  1,4-dinitropiperazine  (34),  1,3-dinitroimida- 
zolidine  (33)  and  RDX  (4)  are  17.1, 18.5  and  19.4,  respectively,  which  show  a  correlation 
coefficient  of  1.00  with  shock  sensitivity  [Murray,  Lane,  Brinck  and  Politzer,  1992]. 
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V.  Strained  Cayelike  Molecules 

Strained  moljn:ules  have  been  a  topic  of  both  fundamental  and  ^lied  interest  in  chemistry 
for  more  than  a  century  [Baeyer,  1885;  Greenberg  and  Liebman,  1978].  The  Commission  on 
Physical  Organic  Chemistry  defines  the  term  strain  as  follows:  "Strain  is  present  in  a  molecular 
entity  or  transition  structure  if  its  energy  is  increased  because  of  unfavorable  bond  lengths,  bond 
angles,  or  dihedral  angles  relative  to  a  standard."  [Doyle,  1992]  In  the  area  of  energetic  materials 
research,  strained  frameworks  are  useful  starting  points  in  designing  new  systems  because  of  their 
higher  energy  content,  which  is  subsequently  released  upon  decomposition. 

Three-,  four-  and  five-membered  rings  (42  -  44)  fused  together  in  a  cagelike  fashion  are 
often  classified  as  strained  "cagelike"  molecules.  Among  hydrocarbons,  some  examples  are 
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tetrahedrane  (45),  triprismane  (46),  cubane  (47)  and  homocubane  (48).  A  key  characteristic  of 
42, 43  and  45 -48  is  the  small  endocyclic  bond  angles  imposed  by  the  three-  and  four-membered 
rings;  for  tetracoordinate  carbon,  normal  bond  angles  are  in  the  neighborhood  of  109°.  This  bond 
angle  distortion  is  a  destabilizing  factor  contributing  to  the  overall  strain  inherent  in  these  systems 
[Greenberg  and  Liebman,  1978]. 

We  have  earlier  developed  an  index  which  provides  an  effective  means  for  characterizing 
individual  strained  bonds  [Politzer,  Abrahmsen,  Sjoberg  and  Laurence,  1983];  it  is  based  on  the 
"bond  path"  concept  [Bader,  Tang,  Tal  and  Biegler-Kdnig,  1982a  and  1982b].  A  bond  path  is  the 
path  between  two  bonded  nuclei  that  goes  through  the  saddle  point  in  the  intemuclear  electronic 
density  and  follows  the  maximum  density  gradient  to  each  nucleus;  it  can  be  visualized  as  the  ridge 
of  highest  electronic  density  linking  the  two  nuclei.  For  many  bonds,  such  as  those  between  the 
carbons  in  propane  and  cyclohexane,  the  bond  path  is  essentially  identical  with  the  intemuclear 
axis;  in  other  cases,  such  as  molecules  with  inherent  bond  angle  strain,  there  can  be  a  significant 
difference  between  the  two  [Politzer,  Abrahmsen,  Sjoberg,  and  Laurence,  1983;  Politzer,  Bar- 
Adon  and  Miller,  1987;  Politzer,  Jayasuriya  and  2^es,  1985;  Politzer,  Kirschenheuter  and  Alster, 


1987;  Politzcr  and  Murray,  1988;  Politzer  and  Seminario,  1989].  Our  "bond  deviation  index"  was 
introduced  as  a  numerical  measure  of  the  degree  to  which  a  given  bond  path  differs  from  a  defined 
reference  path;  it  can  be  viewed  as  reflecting  bond  angle  strain.  A  thorough  discussion  of  this 
index  is  given  elsewhere  [Politzer  and  Murray,  1988]. 

In  Table  5  are  listed  our  bond  deviation  indices  for  a  variety  of  strained  hydrocarbons.  It 
can  be  seen  that  the  C-C  bonds  of  three-membered  rings  have  the  highest  values,  even  when  the 
latter  is  fused  to  a  larger  ring.  Also  included  in  Table  5  are  calculated  isodesmic  reaction  energies 
(computed  at  the  HF/3-21G  level)  for  some  of  the  molecules;  these  values  can  be  interpreted  as 
reflecting  relative  degrees  of  overall  molecular  strain.  The  AEjsodesmic  values  of  the  cagelike 
frameworks  45  -  47  are  considerably  greater  in  magnitude  than  their  respective  building  blocks 
cyclopropane  and  cyclobutane  (as  might  be  anticipated).  However  the  individual  bonds  within  45 
-  47  can  be  characterized  as  either  "cyclopropane-like"  or  "cyclobutane-like"  by  inspection  of  their 
bond  deviation  indices. 

Despite  the  difference  in  the  bond  angle  strain  associated  with  three-  and  four-  membered 
rings,  as  is  reflected  in  their  bond  deviation  indices,  cyclopropane  and  cyclobutane  are  reported  to 
have  very  similar  strain  energies  [Greenberg  and  Uebman,  1978].  The  low,  seemingly  anomalous 
strain  energy  of  the  former  is  attributed  to  the  stabilizing  effect  of  "a-aromadcity"  [Dewar,  1979}. 
This  is  important  to  keep  in  mind  in  designing  potential  energetic  systems.  Cyclopropane  itself 
does  not  contain  a  significant  amount  of  strain  energy;  however  modifications  of  this  ring  that 
disrupt  the  a-aromaticity,  such  as  unsaturation  or  the  formation  of  cagelike  systems,  can  result  in  a 
high  strain  energy. 

Another  interesting  feature  of  strained  hydrocarbons  is  that  regions  of  negative  electrostatic 
potential  ate  associated  with  their  C-C  bonds  [Politzer,  Domelsmith,  Sjoberg  and  Alster,  1982; 
Politzer,  Domelsmith  and  Abrahmsen,  1984;  Politzer  and  Murray,  1988].  (This  is  not  observed 
for  the  C-C  bonds  of  unstrained  hydrocarbons.)  The  most  negative  potentials  are  near  the  bond 
midpoints,  and  their  presence  indicates  that  these  bonds  themselves  can  serve  as  initial  sites  for 
electrophilic  attack.  The  olefin-like  properties  of  cyclopropane  [Cromwell  and  Graff,  1952]  and 
bicyclobutane  [Wiberg,  1965]  and  the  tendency  of  the  former  for  edge  hydrogen  bonding  [Legon, 
Aldrich  and  Flygare,  1982]  can  be  interpreted  as  reflecting  this  feature,  as  can  the  observed  cation- 
catalyzed  rearrangements  of  cubane  and  homocubane  [Cassar,  Eaton  and  Halpem,  1970;  Paquette, 
Boggs  and  Ward,  1975],  in  which  the  inital  step  is  an  interaction  between  the  positive  ion  and  a 
C-C  bond. 

It  is  perhaps  tempting  to  view  the  existence  of  bent  bonds  in  strained  molecules  as  the 
reason  for  the  negative  bond  potentials.  However  we  have  found  that  the  substitution  of  electron- 
withdrawing  groups,  such  as  NO2  and  Q,  eliminates  the  negative  C-C  bond  potentials  in  strained 
hydrocarbons  but  does  not  significantly  affect  the  degrees  of  curvature  of  the  bond  paths  (or  their 
bond  deviation  indices)  [Politzer  and  Murray,  1988].  On  the  other  hand,  density  difference  plots 
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do  show  that  electronic  charge  has  been  depleted  in  the  intemuclear  regions  of  these  derivatives. 
Thus,  two  key  points  that  have  emerged  from  these  findings  are  that  (a)  bent  bonds  do  not 
necessarily  produce  negative  bond  potentials,  and  (b)  the  degree  of  bending  in  a  strained  bond  is 
not  directly  related  to  tli^  quantity  of  electronic  charge  in  the  bond  region. 

We  have  also  shown  that  the  introduction  of  oxygen  or  nitrogen  into  strained  rings 
diminishes  the  bond  deviation  index  of  the  C-X  bond  (where  X  is  O  or  N)  relative  to  the  parent 
hydrocarbon  [Politzer,  Abrahmsen,  Sjoberg  and  Lairence,  1983;  Politzer,  Bar-Adon  and  Miller, 
1987;  Politzer  and  Seminario,  1989],  suggesting  a  slight  decrease  in  bond  angle  strain.  In 
addition,  endocyclic  nitrogens  or  oxygens  effectively  eliminate  any  negative  potentials  associated 
with  C-C  bonds;  instead  characterutic  strongly  negative  regions  are  associated  with  the  nitrogens 
and  oxygens,  indicating  these  to  be  favorable  sites  for  electrophilic  attack  [Politzer  and  Murray, 
1991].  The  slightly  stabilizing  efrect  of  the  introduction  of  oxygen  or  nitrogen  in  strained  ring 
systems  is  also  reflected  in  the  molecule  as  a  whole,  as  can  be  seen  hrom  calculated  isodesmic 
reaction  energies  [Alkorta,  Elguero,  Rozas  and  Balaban,  1990;  Murray,  Seminario,  Lane  and 
Politzer,  1990;  Politzer  and  Murray,  1992];  for  example,  the  calculated  HF/3-21G  AEisodesmic  of 
azacubane  is  121  kcal/mole,  compared  to  127  kcal/mole  for  cubane  (47).  Thus  the  presence  of  aza 
nitrogens  in  cagelike  systems  has  been  suggested  as  being  desirable  [Murray,  Seminario  and 
Politzer,  1991],  in  that  they  may  mitigate  excessive  strain.  For  although  strain  is  viewed  as 
advantageous  for  energetic  systems,  a  high  degree  of  it  can  lead  to  practical  problems  in  the 
synthesis  and  isolation  of  a  proposed  compound. 

In  the  past  decade,  a  variety  of  strained  cagelike  compounds  have  been  proposed  as  target 
energetic  molecules.  These  have  included,  for  example,  nitro  derivatives  of  cubane  and  other  cage¬ 
like  frameworks.  To  further  the  synthesis  of  such  compounds,  we  have  investigated  computation¬ 
ally  the  bond-weakening  tendency  associated  with  NO2  and  NH2  being  bonded  to  adjacent  tertiary 
carbons  in  strained  cagelike  hydrocarbons,  which  has  hampered  efforts  to  synthesize  1,2-dinitro- 
cubane  from  l-amino-2-nitrocubane  [Griffin  and  Marchand,  1989;  Marchand,  1989]. 

When  NO2  and  NH2  are  substituents  on  adjacent  tertiary  carbons  in  some  strained 
molecules,  e.g.  l-nitro-3-atninobicyclobutane  (49)  and  l-amino-2-nitrocubane  (50),  there  is  a 
marked  weakness  of  the  intervening  C-C  bond  (i.e.  C1-C3  in  49  and  C1-C2  in  40),  as  measured 
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by  means  of  the  calculated  bond  order  [Politzer,  Klrschenheuter  and  Alster,  1987;  Murray, 
Seminario  and  Politzer,  1991].  This  phenomenon  is  well  known  to  organic  chemists  and  has  been 
attributed  to  a  "push-pull"  mechanism,  which  in  the  limiting  case  would  lead  to  a  rupture  of  the 
C-C  bond,  as  illustrated  in  eq.  (20): 
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As  an  extreme  example  of  this,  we  have  found  that  the  conformer  of  aminonitrotetrahedrane  in 
which  the  O-N-0  and  N-C1-C2  planes  are  initially  set  to  be  perpendicular  (51)  does  not  optimize  to 
a  local  energy  minimum  but  rather  rearranges  preferentially  to  l-amino-3-nitrocyclobutadiene  (52), 
as  shown  in  eq.  (21)  [Grodzicki,  Seminario  and  Politzer,  1990]. 
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In  more  recent  computational  analyses,  we  have  found  that  the  push-pull  mechanism 
depicted  by  eq.  (20)  is  actually  the  consequence  of  two  separate  effects,  associated  with  the  NH2 
and  NC)2  groups  individually,  which  can  in  some  instances  reinforce  each  other  [Murray, 
Seminario  and  Politzer,  1991a,  1991b;  Murray,  Concha,  Seminario  and  Politzer,  1991].  In 
studies  of  various  aminocubanes,  aminoazacubanes,  aminotriprismanes  and  aminoazatriprismanes, 
we  found  a  consistent  direction-specific  bond  weakening  in  one  (and  only  one)  of  the  C-C  or  C-N 
bonds  adjacent  to  the  site  of  NH2  substitution.  The  weakened  C-X  bond  (where  X=C  or  N)  is 
invariably  coplanar  with  the  C-NH2  bond  and  the  position  of  the  most  negative  electrostatic 
potential  (Voun)  associated  with  the  amine  nitrogen  lone  pair.  This  tendency  for  a  co-planar 
alignment  of  the  lone  pair  and  a  C-X  bond  is  a  manifestation  of  the  anomeric  effect  [Dewar,  1984] 
and  reflects  the  possibility  of  the  delocalization  of  the  lone  pair  as  shown  in  eq.  (22): 
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This  delocalization,  which  weakens  the  C-X  bond,  can  occur  most  effectively  when  the  latter  is 
coplanar  with  the  lone  pair.  We  have  found  that  this  effect  is  stronger  when  X=N  (rather  than 
X=C),  consistent  with  the  anomeiic  effect  being  most  commonly  invoked  for  molecules  in  which 
X  is  an  electronegative  atom,  e.g.  N,0  or  F. 

The  situation  involving  the  NO2  group  alone  is  less  straightforward  than  that  of  the  NH2, 
because  it  can  strengthen  some  bonds  and  weaken  others,  depending  upon  its  conformation 
[Grodzicki,  Seminario  and  Politzer,  1990;  Murray,  Seminario  and  Politzer,  1991a;  Murrray, 
Concha,  Seminario  and  Politzer,  1991].  For  example,  in  smdying  nitrotetrahedrane  (53)  and 
nitrocubane  (54)  we  found  that  when  the  0-N*0  plane  is  perpendicular  to  an  N-C-C  plane,  then 
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that  C-C  bond  is  weakened  slightly;  however  C-C  bonds  that  are  coplanar  (or  nearly  so)  are 
slightly  strengthened.  It  was  shown  through  an  orbital  correlation  analysis  of  53  that  the 
weakening  can  be  attributed  to  an  interaction  between  a  7t-type  NO2  oibital  and  the  C-C  bonding 
orbital  [Grodzicki,  Seminario  and  Politzer,  1990]. 

Our  results  for  nitrotriprismane,  55,  do  not  follow  the  same  pattern  as  was  found  in  53 
and  54  [Murray,  Concha,  Seminario  and  Politzer,  1991].  In  nitrotriprismane  we  found  computa¬ 
tionally  that  only  one  orientation  of  the  NO2  corresponds  to  a  local  energy  minimum;  this  is  when 
the  NO2  group  is  coplanar  with  the  C-C  bond  joining  the  three-membered  tings.  This  is  accom¬ 
panied  by  a  weakening  of  the  adjacent  three-membered-ring  C-C  bonds.  This  was  rationalized  by 
analysis  of  the  electronic  structure  of  triprismane. 

Thus  vicinal  NH2,  NO2  groups  will  weaken  the  C-C  bond  joining  them  when  both  the 
NH2  and  NO2  bond-weakening  tendencies  affect  this  particular  bond  and  reinforce  each  other,  e.g. 
in  49  and  50;  this  can  be  described  by  the  push-pull  mechanism  depicted  in  eq.  (20).  However 
neighboring  NO2  and  NH2  groups  may  also  act  upon  different  bonds,  e.g.  in  the  cases  of  56  and 
57;  in  these,  the  NO2  weakens  the  adjacent  C-C  bonds  in  the  three-membered  ring,  while  the  NH2 
acts  on  the  bond  in  the  four-membered  ring  [Murray,  Concha,  Seminario  and  Politzer,  1991]. 
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In  the  course  of  our  studies  of  nitro  derivatives  of  aza  cagelike  systems,  a  finding  of 
interest  for  the  design  of  new  energetic  materials  has  emerged.  This  is  the  observation  that  the  C-C 
bond  weakening  that  we  have  found  to  be  associated  with  NC)2  in  these  systems  does  not  occur  to 
any  significant  degree  when  there  is  an  adjacent  aza  nitrogen  [Murray,  Seminario  and  Pohtzer, 
1991a;  Murray,  Concha,  Seminario  and  Politzer,  1991].  Taken  in  conjunction  with  the  overall 
stabilizing  influence  of  aza  nitrogens  (through  diminishing  the  molecular  strain  energy)  [Murray, 
Seminario,  Lane  and  Politzer,  1990],  this  further  emphasizes  the  advantages  of  aza  nitrogens  from 

NO2 


the  standpoints  of  synthesis  and  stability.  Thus,  we  have  suggested  that  target  systems  such  as  the 
tetranitrotetraazacubane  58  are  plausible  from  a  synthetic  standpoint  and  should  be  further  investi¬ 
gated  as  potential  high-energy  molecules  [Murray,  Seminario  and  Politzer,  1991a].  (58  could  in 
principle  be  made  firom  amine  precursors,  with  no  possibility  of  vicinal  NH2,N02  combinations.) 
This  suggestion  is  generally  consistent  with  comments  made  by  Engelke  and  Stine  'manab  initio 
SCF-MO  study  of  octaazacubane  (Ng)  [Engelke  and  Stine, 1990];  they  stated  that  "the  known  high 
stability  of  (CH)g  and  the  predicted  high  energy  of  Ng  cubane  suggest  that  there  are  probably 
intermediate  mixed  carbon-nitrogen  cubane-like  structures  that  would  be  stable  high-energy 
materials."  This  comment  was  applied  also  to  nitrated  azacubanes. 

VI.  High-Nitrogen-Content  Molecules 

There  is  a  growing  interest  in  high-nitrogen-content  molecules  as  potential  energetic 
materials  [Benson,  1984],  prompting  a  variety  of  theoretical  studies  [e.g.  Mayers,  Salter,  Bartlett, 
Salter,  Hess  and  Schaad,  1988;  Engelke,  1989, 1990;  Engelke  and  Stine,  1990;  Politzer  and 
Seminario,  1990;  Thomas,  Quelch  and  Schaefer,  1991;  Seminario,  Concha  and  Politzer,  1992; 
Politzer,  Murray,  Seminario  and  Miller,  1992].  In  this  section  we  will  discuss  calculations  that 
have  been  carried  out  for  molecules  in  this  category  which  contain  three  or  more  linked  nitrogens 
(nitrogen  catenation). 

Nitrogen  catenation  is  often  associated  with  instability  and/or  sensitivity  [Benson,  1984]. 
For  example,  the  powerful  explosive  l,3,5-triazido-2,4,6-trinitrobenzene  (59)  undergoes  a  slow 
conversion  to  hexanitrosobenzene  (60)  and  N2  [eq.  (23)], 


25 


eventually  losing  its  initiating  power  [Meyer,  1987].  On  the  other  hand,  the  tetrazene  61  has  been 
found  to  be  quite  insensitive  to  impact  [Benson,  1984]  and  the  tetraazapentalene  62  is  known  to 
have  a  high  level  of  thermal  stability  [Meyer,  1987].  These  findings  raise  the  question  of  what 


factors  are  important  in  determining  the  stabilities  of  molecules  with  nitrogen  catenation.  Several 
studies  addressing  this  point  will  be  discussed. 

In  an  earlier  computational  analysis  of  some  tetraaz^ntalenes  [Politzer  and  Seminario, 
1990],  we  found  that  the  parent  molecule  63  is  significantly  more  stable  than  the  isomer  64.  Our 
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computed  structure  for  63  showed  that  there  is  considerable  delocalization  of  charge,  as  can  be 
depicted  by  resonance  structures  63a  -  63e.  In  contrast,  64  has  localized  N=N  and  C=C  double 
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bonds.  The  former  may  be  a  destabilizing  feature,  since  it  may  promote  decomposition  through 
elimination  of  N2. 

The  in^rtance  of  this  factor  is  brought  out  very  strikingly  by  the  pairs  of  isomeric 
triazoles  65, 66  and  67, 68;  their  measured  sensitivities,  given  as  impact  drop  heights,  are  shown 
in  parentheses  [Storm,  Ryan,  Ritchie,  Hall  and  Bachrach,  1989].  65  and  67  are  clearly  more 


67  (9  cm)  68  (67  cm) 


sensitive  to  impact  than  their  isomers  66  and  68.  The  crystallographically-determined  structures 
of  65  -  68  show  that  the  primary  structural  difference  within  each  pair  of  isomers  is  that  the  Nb- 
Nc  bonds  in  the  triazole  rings  of  65  and  67  are  about  0.03  A  shorter  than  the  other  N-N  bonds. 

(In  65  and  67,  Na  is  the  nitrogen  linked  to  the  phenyl  portion.)  It  was  suggested  that  this 
facilitates  the  decomposition  of  65  and  67  through  the  loss  of  N2,  accounting  for  the  greater 
sensitivities  of  these  two  isomers  [Storm,  Ryan,  Ritchie,  Hall  and  Bachrach,  1989]. 

The  evidence  presented  suggests  that  nitrogen  catenation  may  be  compatible  with  stable 
systems  provided  that  there  is  considerable  conjugation  between  the  nitrogens;  one  indication  of 
this  is  similarity  of  the  N-N  bond  lengths.  The  presence  of  some  shorter,  stronger  N-N  bonds  is  a 
warning  of  possible  instability/sensitivity  due  to  f^ile  decomposition  through  N2  elimination.  (An 
extreme  example  is  provided  by  the  azides,  R-N=N=N,  many  of  which  are  highly  explosive,  and 
decompose  through  release  of  N2;  the  outer  N-N  bond  lengths  are  typically  about  1.12  A  [Allen, 
Keimard,  Watson,  Brammer,  Oipen  and  Taylor,  1987],  nearly  the  same  as  free  N2, 1.098  A.) 
This  structural  criterion  can  also  be  applied  to  a  recent  computational  study  of  the  unknown 
unsubstituted  tetrazines  69  and  70  [Thomas,  Quelch  and  Schaefer,  1991}.  In  69,  the  ab  initio 
CISD/DZP  computed  N-N  bond  lengths  are  1.330, 1.297  and  1.330  A;  in  70  they  are  both 
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1.31 1  A.  70  is  found  to  be  the  more  stable  of  the  two  isomers.  The  only  known  tetrazine  is  71, 
in  which  the  CISD/DZP  N-N  distances  are  1.310  A. 

In  the  case  of  hypothetical  molecules  that  have  not  yet  been  observed  or  isolated,  a  key 
question  is  whether  they  can  exist  even  in  principle.  This  can  be  addressed  computationally  by 
carrying  out  a  vibrational  frequency  analysis  for  the  structure  of  interest;  if  this  yields  no  imaginary 
frequencies,  then  the  stmcture  corresponds  to  a  local  minimum  on  the  energy  surface  [Hehre, 
Radom,  Schleyer  and  Pople,  1986].  For  example,  the  unknown  tetrazines  69  and  70  have  both 
been  shown  computationally  to  be  local  minima  on  the  C2N4H2  energy  surface  [Thomas,  Quelch 
and  Schaefer,  1991].  The  same  is  true  of  pentazine  (72)  [Politzer,  Murray,  Seminario  and  Miller, 
1992],  and  the  Ng  analogue  of  benzene  (73),  along  with  four  of  its  cyclic  isomers  74  -  77 
[Engelke,  1989],  none  of  which  have  as  yet  been  confirmed  experimentally.  In  the  case  of  72,  the 
computational  analysis  was  at  the  post-Hartree-Fock  MP2/6-31G*  level,  but  for  73  -  77,  the 
conclusions  were  drawn  from  HF/4-31G*  calculations.  73  is  predicted  to  be  the  most  stable  of  the 
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cyclic  Ne  isomers  [Engelke,  1989],  which  is  consistent  with  its  having  six  equivalent  N-N  bonds 
and  aromatic  character,  74, 75  and  77  each  have  one  or  two  localized  N-N  double  bonds,  while 
hexaazatriprismane  (76)  has  two  types  of  N-N  single  bonds. 

It  has  been  demonstrated  by  us  [Seminario,  Concha  and  Politzer,  1991]  and  by  others 
[Clabo  and  Schaefer,  1987;  Coffin  and  Pulay,  1991]  that  reliable  assessment  of  the  structures  and 
stabilities  of  molecules  containing  a  large  proportion  of  strongly  electronegative  atoms  (N,  O,  F, 
etc.)  requires  that  electronic  correlation  be  taken  into  account,  especially  if  these  atoms  are  linked, 
as  in  H2O2  [Politzer  and  Murray,  1992].  For  example,  in  a  study  of  fiiroxan  (78:  1,2,5- 
oxadiazole-2-oxide)  and  some  of  its  isomers  (e.g.  79  -  82),  we  showed  that  the  structural  data 
and  predicted  stabilities  were  markedly  dependent  upon  the  computational  method  and  level  used 
for  the  geometry  optimizations.  For  instance,  MP2/6-31G*  calculations  revealed  81  and  82  to  be 
unstable,  although  HF/3-21G  and  HF/6-31G*  had  shown  local  minima. 
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Pentazine,  72,  is  predicted  to  be  a  stable  system  at  both  the  HF/3-21G  and  MP2/6-3 IG* 
levels;  the  computed  geometries  are  very  similar,  with  the  C-N  and  N-N  bond  distances  differing 
by  0.014  and  less  than  0.009 A,  respectively  [Politzer,  Murray,  Seminario  and  Miller,  1992].  On 
the  other  hand,  the  N-N  bond  lengths  calculated  by  post-Hartree-Fock  methods  for  triaziridine 
(83)  are  0.054  A  longer  than  those  obtained  at  the  HF/6-31G*  level  [Magers,  Salter,  Bartlett, 
Salter,  Hess  and  Schaad,  1988].  It  may  be  that  Hartree-Fock  methods  are  reliable  for  aromatic  aza 
systems,  e.g.  pyridine,  the  tetrazines,  pentazine  and  N6  benzene,  but  they  certainly  should  be 
applied  with  caution  to  high-nitrogen  systems  in  general. 


Our  shock  sensitivity  relationships  for  nitramines  and  nitroaliphatics  suggest  that  the 
strengths  and  numbers  of  X-NO2  linkages  in  relation  to  overall  molecular  size  are  important  in 
determining  sensitivity  [Politzer,  Murray,  Lane,  Sjoberg  and  Adolph,  1991].  Our  observation,  for 
a  few  nitramines,  that  sensitivity  varies  directly  with  local  polarity  (as  measured  by  IT)  is  consistent 
with  this  relationship  [Murray,  Lane,  Brinck  and  Politzer,  1992].  It  is  intuitively  reasonable  to 
anticipate  that  local  polarity  should  increase  as  the  number  of  X-NO2  linkages  per  unit  volume 
increases. 

In  the  area  of  nitroaromatics,  correlations  between  impact/shock  sensitivity  and  total 
electron-donating  ability  of  the  non-nitro  substituents  [Owens,  1984],  or  Vmid  or  Vmidjnax 
[Owens,  Jayasuriya,  Abrahmsen  and  Politzer,  1985;  Murray,  Lane,  Politzer  and  Bolduc,  1990] 
suggest  that  the  strengthening  of  C-NO2  linkages  leads  to  decreased  sensitivity.  A  direct 
relationship  between  nitroaromatic  impact  sensitivity  and  local  polarity  [Murray,  Lane  and  Politzer, 
1992],  is  again  not  inconsistent  with  this  viewpoint,  in  that  11  decreases  in  response  to  the  addition 
of  electron-donating  groups  and  increases  with  the  addition  of  electron-withdrawing  groups.  The 
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implications  of  correlations  between  local  polarity  and  impact/shock  sensitivity  are  being  further 
investigated. 

Molecules  which  do  not  fit  the  general  types  of  relationships  that  have  been  found  to  exist 
between  impact/shock  sensitivity  and  molecular  and/or  structural  properties  suggest  starting  points 
for  further  investigations  aimed  at  understanding  critical  molecular  features  which  may  be 
determinants  of  sensitivity  in  certain  specific  cases.  For  example,  we  have  suggested  that  the 
anomalously  high  sensitivities  of  hydroxyaromatics  may  be  due  to  formation  of  small  quantities  of 
unstable  nitronic  acid  tautomers  [Politzer,  Seminario  and  Bolduc,  1989].  Table  4  shows  that  the 
surface  electrostatic  potential  of  the  OH  hydrogen  in  the  monohydroxy  derivative  picric  acid  13  is 
considerably  more  positive  than  the  areas  above  and  below  the  aromatic  ring.  This  is  unusual  for 
trinitroaromatics,  and  supports  our  earlier  specul^on  that  the  facile  loss  of  a  phenolic  hydrogen 
may  in  part  account  for  the  high  sensitivities  of  hydroxynitroaromatics. 

The  l-picryl- 1,2,3-triazoles  65  and  67  are  another  example  of  anomalously  high 
sensitivity  being  attributed  to  a  molecular  feature  other  than  X-NC)2  linkages,  in  particular  the 
presence  of  a  relatively  short  N-N  bond  which  may  favor  the  facile  elimination  of  N2  [Storm, 
Ryan,  Ritchie,  Hall  and  Bachrach,  1989].  This  may  also  be  an  important  factor  in  determining  the 
sensitivities  of  the  isomers  4-nitro- 1,2,3-triazole  (84)  and  3-nitro- 1,2,4-triazole  (85);  the  impact 
drop  heights  of  these  are  25  and  >320  cm,  respectively.  We  have  found  that  the  former,  more 
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impact  sensitive  compound  has  calculated  HF/3-21G  N-N  bond  lengths  of  1.287  and  1.398  A, 
respectively,  while  the  N-N  distance  in  85  is  1.420  A.  The  relatively  short  N-N  bond  length  in  84 
is  suggestive  of  a  localized  double  bond,  favoring  the  elimination  of  N2. 

In  general  -litriurines  are  more  sensitive  to  impact/shock  than  are  nitroaromatics;  e.g.  see 
Table  1.  In  this  contcru  it  may  be  relevant  to  point  out  that  the  experimentally-determined  C-NO2 
rotational  energy  barriers  for  nitromethane,  nitrobenzene  and  nitroethylene  range  from  0-5 
kcal/mole  (increasing  in  that  order)  [H0g,  Nygaard  and  Sorensen,  1971;  Meyer,  Gammeter, 

Smith,  Kiihne,  Ndsberger  and  Giinthard,  1973;  Conel,  Larson  and  Pedersen,  1980],  while  that  of 
the  N-NO2  bond  in  dimethylnitramine  is  greater  than  9  kcal/mole  [Kintzinger,  Lehn  and  Williams, 
1969],  as  are  also  the  theoretically  predicted  N-NO2  rotational  barriers  of  a  series  of  alicyclic 
nitramines  [Habibollahzadeh,  Murray,  Redfera  and  Politzer,  1991;  Habibollahzadeh,  Murray, 
Grice  and  Politzer,  1992].  It  was  proposed  by  Kamlet  and  Adolph  [1981]  that  X-NC)2  rotation 
provides  a  means  for  dissipating  thermal  energy  that  otherwise  might  be  used  to  initiate 
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decomposition.  Thus  it  can  be  speculated  that  the  lower  rotational  energy  barriers  of  C-NO2  bonds 
in  nitroaromatics  relative  to  the  N-NC)2  in  nitramines  may  contribute  to  the  generally  higher 
sensitivities  of  the  latter.  The  higher  barrier  to  rotation  of  the  N-NC)2  linkage  may  impede  the 
dissipation  of  thermal  energy  and  lead  to  generally  greater  sensitivity. 

Our  discussion  of  serrsitivity  and  its  relationship  to  molecular  properties  has  focused  upon 
general  correlations  that  exist  within  separate  classes  of  compounds,  e.g.  nitroaromatics, 
nitramines  and  nitroaliphatics,  and  upon  exceptions  to  these  relationships  which  suggest  alternative 
mechanisms  to  X-NO2  bond  cleavage  as  being  operative  in  specific  instances.  Thus,  although  we 
believe  it  to  be  unlikely  that  an  unequivocal  universal  relationship  between  sensitivity  and 
molecular  parameters  can  be  established,  considerable  progress  has  indeed  been  made  in 
understanding  factors  that  influence  sensitivity  within  separate  classes  of  energetic  systems  as  well 
as  in  certain  exceptional  cases.  In  the  area  of  high-nitrogen-content  molecules,  it  spears  that 
reliable  structural  analyses  are  particularly  important  in  assessing  stability  and/or  sensitivity;  in  a 
number  of  cases,  the  presence  of  localized  N-N  bonds  appears  to  be  a  factor  favoring 
instability/high  sensitivity. 
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Table  1.  Measured  impact  and  shock  sensitivities  of  some  common  explosives.^ 


Molecule 


Impact  drop  height^ 
(cm) 


(kbar) 


O2N— O-HjC  CH2-O-NO2 
02N-0-H2C^  ^CH2-0-N02 
1,  PETN 


CH2-O— NO2 
CH-0— NO2 
CHz-0— NO2 
2,  nitroglycerin 


r 'I 

N  N 

O2N'  'NO2 


10.97 


4,  RDX 


O2N  ^  NO2 

'n  n" 

<  > 

N  N 

02N^  ^  'NO2 


10.81 


7,  HMX 


HaC^  NO2 
N 


10.64 


5,  Tetryl 


(continued) 


Table  1.  Measured  impact  and  shock  sensitivities  of  some  common  explosives  (continiied).^ 


Molecule 


Impact  drop  height^’ 
(cm) 


(kbar) 


3,  TNB 


8.  TNT 


100 


160 


14.96 


17.72 


9,  MATE 


177 


27.91 


320 


46.20 


_ 10,  TATB _ 

^Data  taken  from  Storm,  Stine  and  Kramer,  1990. 

•’The  impact  sensitivity  as  impact  drop  height  as  measured  at  the  Los  Alamos  National  Laboratory 
and/or  Naval  Surface  Warfare  Center  using  the  Bruceton  Method,  type  12  tools,  2.5  kg  weight,  40 
mg  sample,  S/0  sandpaper,  25  trials.  Impact  sensitivity  increases  as  impact  drop  height  decreases. 

'’The  shock  sensitivity  as  measured  by  the  Naval  Surface  Warfare  Center  small  scale  gap  test;  P90  is 
the  pressure  to  90%  of  its  theoretical  maximum  density.  Shock  sensitivity  increases  as  P90 
increases. 
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Table  3.  Sonae  calculated  and  experimental  moleoiiar  ^gppemes.^ 


Mdecule 

n 

(kcal/mole) 

(Debye) 

e 

C-C6H12 

2.2 

0.0 

2.02 

C6H5CH3 

4.6 

0.36 

2.38 

CQ4 

4.8 

0.0 

2.24 

C6H6 

4.9 

0.0 

2.28 

CeHscr  H3 

7.2 

3.02 

17.39 

CeHsa 

8.0 

1.69 

5.71 

C6H5OH 

8.6 

1.45 

9.78 

CH3COCH3 

9.4 

2.88 

20.7 

C2H5OH 

10.1 

1.69 

24.3 

C6H5NO2 

12.3 

4.22 

34.82 

CH3C00H 

12.9 

1.74 

6.15 

P-C6H4(N02)2 

16.5 

0.00 

— 

CH3CN 

17.1 

3.92 

36.2 

(CH3)2S0 

19.1 

3.96 

46.60 

CH3NO2 

19.9 

3.46 

35.9 

H2Q 

21,6 

1.85 

78.54 

^Tbese  calculated  11  values  and  experimentally  determined  dipole  moments  \i  and 
constants  e  are  taken  from  Biinck  et  al  [1992]. 
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Table  4.  Calculated  and  experimentally-deteronned  properties  fey  some  trinitroaromatics. 


X 

O^N^L  NO2 

N02 

X,Y.Z 

iopict 
drop  height^ 

(cm) 

Vjnidjnax^ 

(kcal/mole) 

nc 

(kcal/mole) 

Ring  center 
Vs.max*^ 
(kcal/mole) 

H,  NO2,  H 

— 

142.4 

21.4 

52 

1,2,4, 5-tetranitro 

27 

141.8 

21.4 

50 

NH2,  NO2.  H 

41 

150.6 

21.2 

46 

NH2:  42 

Cl,  H,  H 

79 

146.8 

19.1 

45 

OH,  H,  H 

87 

127.4 

21.0 

OH:  58 

39 

H,H,  H 

100 

130.5 

19.5 

39 

COOH,  H,  H 

109 

128.6 

19.5 

COOH:  43 

41 

CH3.H,H 

160 

131.1 

18.2 

37 

NH2,H,H 

177 

109.2 

18.6 

36 

NH2:  32 

OCH3,  H,  H 

192 

117.3 

18.2 

40 

OCH3,  OCH3,  H 

251 

116.1 

17.6 

38 

NH2,  NH2,  H 

320 

91.6 

17.1 

33 

NH2,  NH2,  NH2 

>320 

64.6 

16.1 

31 

arfhe  impact  sensitivity  as  impact  drop  height  as  measured  at  the  Los  Alamos  National 
Laboratoiy  and/or  Naval  Si^ace  Waifate  Center  using  the  Bruceton  Metho^  type  12 
tools,  2.S  kg  weight,  40  mg  sample,  5/0  sandpaper,  2S  trials.  Impact  sensidvi^  increases 
as  impact  drop  height  decreases. 

I’Data  talfwn  Murray,  Lane,  Politzer  and  Bolduc,  1990;  units  have  been  converted  to 
kcalAxiole. 

^ata  taken  from  Murray,  Lane,  Brinck  and  Politzer,  1992. 

^  The  Vs  jnax  are  located  above  or  below  the  aromatic  ring  unless  otherwise  indicated. 


x:  0.084 
y:  0.080 


0.116 


151.0 


w:  0.023 
x:  0.022 
y:  0.081 
z:  0.078 


x:  0.033 
y:  0.080 


140.3 


0.029 


127.0 


x:  0.030 
y:  0.020 
z:  0.009 

X 

aTaken  from  Politzer  and  Muzcay,  1988. 

^HF/3-21G  results  taken  from  Mmray,  Lane,  Seminario  and  Politzer,  1990,  unless  otherwise 
indicted. 

<^HF/3-21G  results  taken  from  Murray,  Redfem,  Seminario  and  Politzer,  1990. 


Figure  1.  Calculated  HF/STO-SG  electrostatic  potential  on  tibe  molecular  surface  of  benzene. 

Two  ranges  of  V(r)  are  depicted;  these  are,  in  kcalAnole:  white  <  0;  0  <  gray  <  10. 

Figure  2.  Calculated  HF/STO-SG  electrostatic  potential  on  die  molecular  surface  of  nitrobenzene. 
Three  ranges  <kv(r)  are  dqncted;  th^  are,  in  kcalAnole:  white  <  0;  0  <  gray  <  10; 

10  <  blac^ 

Figure  3.  Calculated  HF/STO-SG  electrostatic  potential  on  the  molecule  surface  of  1,3,S- 

trinitrobenzene.  Three  ranges  of  V(r)  are  depicted;  these  are,  in  kcalAnole;  white  <  0; 

0  <  gray  <  3S;  3S  <  black. 


Figure  1.  Calculated  HF/ST0-5G  electrosuuic  potential  on  the  molecular  surface  of  benzene. 

Two  ranges  of  V(r)  are  depicted;  are,  in  kcal/mole:  white  <  0;  0  <  gray  <  10. 


Hgure  2.  Calculated  HF/STO-5G  electrostatic  potential  on  the  molecular  surface  of  nitrobenzene. 
Three  ranges  of  V(r)  ate  depicted;  th^  are,  in  kcal/mole:  white  <  0;  0  <  gray  <  10; 

10  <  black. 


